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Electromagnetic  (EM)  actuators  in  space  applications  are  not  a  new  idea,  but  they  are  most  of 
the  time  associated  to  low  Earth  orbit  missions,  where  the  on-board  magnetic  moment 
interacts  with  the  Earth  magnetic  field.  More  recently  EM  actuators  have  been  studied  in  the 
context  of  formation  flying  as  a  way  to  generate  inter-spacecraft  force  and  torque  in  order  to 
control  the  formation  geometry  (relative  position  and  attitude  of  the  spacecrafts).  With  respect 
to  other  possible  actuators  (like  EEEP  and  cold-gas  thrusters),  EM  actuators  do  not  require 
propellant  (hence  ensuring  longer  lifetime)  and  do  not  cause  contamination  problems,  but  are 
effective  only  in  a  limited  range  of  distances  (due  to  the  magnetic  field  decreasing  as  1/r  ). 

In  this  paper  we  consider  a  setup  where  a  spacecraft  (called  the  hub)  generates  the  magnetic 
field,  while  the  other  spacecrafts  (called  the  flyers)  modify  their  on-board  magnetic  moment 
in  order  to  generate  the  desired  force  and  torque.  We  show  that  if  the  magnetic  field  generated 
by  the  hub  is  constant,  it  is  not  possible  to  generate  any  combination  of  force  and  torque 
(there  are  forbidden  directions,  independently  of  the  magnitude  of  the  desired  force  and 
torque)  and  we  present  a  possible  solution  based  on  a  rotating  field.  The  main  idea  is  to  have  a 
time-varying  magnetic  field  generated  by  the  hub  (chosen  as  a  rotating  magnetic  dipole)  and 
time- varying  magnetic  moments  on  the  flyers.  The  variation  law  on  each  flyer’s  magnetic 
moment  is  computed  in  order  to  obtain  the  desired  force  and  torque  on  average  (over  one 
rotation  period).  With  such  approach,  the  instantaneous  force  and  torque  may  differ  from  the 
desired  ones,  these  discrepancies  being  “filtered  out”  by  the  spacecraft  inertia. 

Einally  we  present  a  preliminary  sizing  of  the  actuators  based  on  classical  and 
superconducting  coils,  for  typical  mission  requirements  :  this  underlines  the  technological 
challenges  to  be  faced  in  order  to  build  the  proposed  kind  of  actuator. 
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1  Magnetic  actuators  basics 


1 . 1  Notations 

A  given  vector  quantity  will  be  denoted  in  bold  (x)  in  the  text  and  with  an  overhead  arrow  in 
the  equations.  Its  module  will  be  plain  (x)  in  the  text  and  in  the  equations  and  the 
corresponding  unitary  vector  will  be  noted  Ux  (with  an  overhead  arrow  in  the  equations).  Dot 
and  cross  products  are  noted  respectively  •  and  x .  The  vacuum  permeability  is  noted  po  and  is 
4pl0'’  H/m. 

1 .2  Dipolar  field  and  interactions 

Considering  a  dipolar  magnetic  moment  M  located  in  the  origin  of  a  given  reference  frame, 
the  magnetic  field  at  the  point  r  is  given  by  : 

B{M,7)  =  ^\  [3(M  •  F)?  -  r  ^  ^  ] 

4n  r  4n  r 

If  a  second  dipolar  magnetic  moment  is  located  at  r,  the  interaction  force  and  torque  are  : 
f  (in, M ,r)  =  mx  B(M ,  r) 

F(m,M,r)  =  grad  (in  -  B(M  ,r))  =  (m-  grad)B(M  ,r) 

Developing  the  expressions  gives  : 

Po  3(M  ■  r)(mxr)  -  (ihxM) 

4n 
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- ^  M  )\ 
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It  is  obvious  that  given  a  value  of  B  it  is  not  always  possible  to  chose  a  value  of  m  that 
generates  any  desired  value  of  T  and/or  F.  When  M  and  m  are  oriented  as  r,  B  is  oriented  as 
r  too,  the  torque  is  null  and  the  module  of  the  force  (repulsive  if  the  directions  of  M  and  m 
are  opposite  each  other,  attractive  if  they  are  equal)  is  : 

_  3po  Mm 

1 .3  Two  space  applications 

Magneto  torque  bars  (MTBs)  are  frequently  used  to  control  the  attitude  of  low  Earth  orbit 
satellites  :  the  on-board  magnetic  moment  is  modified  by  a  control  law  and  interacts  with  the 
Earth’s  magnetic  field.  The  interaction  torque  is  used  to  control  the  attitude  of  the  satellite  (or 


to  dump  the  rotation  speed  of  the  satellite  after  separation  from  the  launcher).  The  magnetic 
field  seen  by  the  satellite  is  time-varying  (due  to  the  relative  motion  between  the  Earth  and 
the  satellite)  and  different  control  laws  take  this  into  account  in  different  ways. 

In  the  NASA  Terrestrial  Planet  Finder  (TPF)  mission,  n  satellites  are  aligned  and  the  entire 
formation  is  rotating  around  the  central  point :  radial  forces  between  satellites  are  necessary  to 
compensate  the  centrifugal  force  and  to  ensure  relative  positioning.  Magnetic  actuators  have 
been  proposed  for  such  task  and  dedicated  control  laws  have  been  studied.  It  has  been  shown 
that  the  magnetic  actuators  could  also  be  used  to  spin  the  formation  :  in  this  case  however,  the 
magnetic  moments  are  not  co-linear  with  respect  to  the  relative  position  vectors,  hence  each 
spacecraft  sees  a  disturbance  torque  that  must  be  compensated  to  avoid  attitude  modifications. 


2  Proposed  solution 
2.1  Concept 

The  proposed  solution  consists  in  a  hub  satellite  generating  a  time-varying  magnetic  field  and  N  flyer 
satellites.  The  magnetic  moment  of  each  flyer  is  computed  in  order  to  produce  the  desired  torque  and 
force  on  average,  knowing  the  variation  law  of  the  magnetic  field  generated  by  the  hub  and  the  relative 
position  of  the  flyer  with  respect  to  the  hub. 

In  this  paper  we  focus  on  one  of  the  simplest  implementation  :  the  hub  is  located  in  the  origin  of  the 
reference  frame  and  carries  a  magnetic  momentum  M(t),  constant  in  module  but  rotating  around  the 
origin  at  constant  rotation  rate  ?  ,  while  a  single  flyer  located  at  r  carries  the  magnetic  moment 
m(t),  varying  in  a  synchronous  way.  In  Cartesian  coordinates  : 


cos((Ot) 

mc^  cos((£)t)  -1-  ms^  sin((Of) 

M(t)  =  M 

sin((oO 

0 

=  M(M,(p,t)  , 

m(t)  = 

mCy  cos((00  -1-  mSy  sin((Of) 
mc^  cos((Ot)  -1-  ms^  sin((Of) 
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The  parameters  of  the  variation  law  of  the  flyer’s  moment  are  mcx/y/z  ans  msx/y/z  and  are  noted 
p  in  the  rest  of  the  document.  Plugging  the  expressions  of  M(t)  and  m(t)  in  the  expressions  of 
B,  T  and  F  we  get : 

5(f,M(0)  =  — ^[3(M(0-r)f-r"M(0l=5(F,M,(D,0 
4%  r^ 

f(m(t),B(t))  =  m(t)xB(t)  =  f(p,r,M,(0,t) 

F(m(t),B(t))  =  grad(m(t)  ■  B(t))  =  F(p,r,M,(0,t) 

The  above  instantaneous  torque  and  force  are  then  analytically  integrated  over  one  period  to 
obtain  the  expression  of  their  average  value  : 

-  1  ^  - 

Ta  =  —^f(p,r,M,(0,t)dt  =  f^(p,r,M) 

-  1  ^  - 

Fa  =  —^F(p,r,M,(0,t)dt  =  F^(p,r,M) 

^  0 


It  turns  out  that  the  above  expressions  of  the  average  torque  and  force  can  be  rewritten  as  : 


where  C(r,M)  is  a  6  by  6  matrix  depending  only  on  r  and  M.  It  is  at  this  point  clear  that  the 
above  relation  can  be  used  to  compute  p  given  an  average  required  torque  and  force  : 


p(0  =  C-'(F(0,M)- 


The  above  relationship  is  only  valid  where  C  is  non-singular  and  one  can  compute  (the 
Cartesian  components  of  r  being  x,  y  and  z)  : 


det(C(F,M)) 


21M\x^  +  y^)z 
64(v"  +  y"  +  z")" 


The  above  expressions  allows  to  deduce  the  singular  configurations,  namely  all  the  points  in 
the  XY  plane  (the  plane  in  which  lies  the  M(t)  vector)  and  all  the  points  of  the  Z  axis  (the  axis 
of  rotation  of  M(t)).  Close  to  these  singular  configurations  there  are  combinations  of  required 
average  force  and  torque  that  needs  an  extremely  large  magnetic  momentum  to  be  obtained. 


2.2  Remarks 

In  presence  of  multiple  flyers  the  situation  does  not  change  too  much,  but  it  is  necessary  to 
avoid  inter-flyer  magnetic  interaction.  Since  the  maximum  interaction  force/couple  depends 
on  the  product  M-m,  it  is  possible  to  put  small  magnetic  moments  on  the  flyers  and  a  big  one 
on  the  hub  (this  assumes  the  flyer-flyer  distance  to  be  comparable  to  the  hub-flyer  one). 

Notice  also  that  the  hub  will  see  the  sum  over  all  the  flyers  of  the  interaction  torque  and  force 
and  it  will  move  and  rotate  accordingly  :  inertial  actuators  (like  thrusters)  are  hence  required 
to  compensate  these  movements,  but  they  can  be  much  more  coarse  since  the  only  objective  is 
to  stay  away  from  singular  configurations.  This  because  the  mission  requirement  is  to 
maintain  relative  attitude  and  positioning  of  the  flyers,  the  hub  being  not  part  of  the  scientific 
instrument  of  the  mission. 

All  the  computations  rely  on  the  exact  knowledge  of  the  magnetic  field  :  the  dipolar 
approximation  used  is  accurate  only  in  far  field,  i.e.  when  r  is  much  bigger  than  the 
characteristic  length  of  the  magnetic  dipole  that  generates  the  field.  Otherwise  the  average 
torque/force  may  be  affected  by  an  error  (that  could  be  compensated  by  the  outer  control 
loop).  Notice  also  that  the  system  requires  the  knowledge  of  hub-flyer  relative  position  and 
attitude  and  a  synchronisation  (the  flyer’s  magnetic  moment  law  must  be  synchronised  to  the 
hub  magnetic  moment  rotation). 


3  Preliminary  sizing 

In  the  following  paragraphs  we  will  show  some  very  preliminary  computations  in  order  to 
highlight  the  technological  challenges  involved  in  the  realisation  of  the  proposed  solution. 


3.1  Overall  sizing 

Supposing  that  the  magnetic  moment  of  the  flyer  is  generated  by  3  orthogonal  moments,  we 
focus  on  one  of  them  (let’s  take  the  x  axis)  and  we  obtain  the  following  relations  between  the 
2-norm  of  the  averaged  desired  torque/force  and  the  x-axis  magnetic  momentum  of  the  flyer  : 


me. 

_ 
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The  C(r,M)  matrix  is  affine  in  M  (C(r,M)  =  M-C(r,l))  :  the  first  quantity  we  want  to  size  is 
the  product  m-M  with  respect  to  the  hub-flyer  distance  and  the  required  force/torque.  As 
stated  before;  when  the  two  magnetic  moments  are  aligned,  the  interaction  force  module  is  : 

_  dp-o  Mm 

This  value  is  an  upper  bound  with  respect  to  what  we  can  obtain  with  our  solution  because  we 
are  dealing  with  non-aligned  moments  (at  least  most  of  the  time)  and  because  we  aim  at 
producing  force  and  torque.  In  fact,  we  can  still  use  the  above  formula  introducing  a 
correction  factor  Cf  : 


3po  Mm  1 
In  cl 

The  correction  factor  depends  on  the  distance  and  on  the  flyer  elevation  angle,  but  for  a  given 
angle  is  pretty  constant  with  respect  to  the  distance  (it  becomes  large  in  the  near  field,  but 
there  our  computations  are  no  longer  valid)  :  we  can  consider  it  to  be  in  the  [4. 0-5. 5]  range  for 
an  elevation  angle  in  [30°-75°],  the  minimum  (and  hence  the  optimum)  being  roughly  located 
at  60°. 


11  mM 


10 


As  an  example,  we  choose  a  flyer  located  at  50m  from  the  hub  and  at  an  elevation  angle  of 
60°.  We  would  like  to  obtain  an  average  torque  and  force  vector  with  a  2-norm  equal  to  1  -6. 

E  2  2 

Our  computations  give  a  value  of  the  m-M  product  of  roughly  1.  8  (A-m  )  and  supposing  the 
M/m  ratio  equal  to  100,  we  obtain  M  =  1  .^5  A-m^  and  m  =  1  .^3  A-m^. 


3.2  Flyer  magnetic  moment  sizing 

The  magnetic  moment  of  the  flyer  is  smaller  and  is  time  varying  :  because  of  that  we  choose 
“classical”  coils  (and  not  superconducting  ones),  with  or  without  a  magnetic  core. 

3.2.1  Flat  coil  without  ferromagnetic  core 

For  a  single  coil  of  radius  R  and  with  a  current  I  obtained  as  the  product  of  the  coil  wire 
section  by  the  current  density  J  (I  =  J-Sc)  the  magnetic  moment  is  : 

M  ={%R^){JSc) 


In  fact,  we  suppose  that  the  coil  wire  with  section  Sc  is  built  with  N  thinner  wires  with  section 
Sw  and  a  filling  ratio  of  1  (optimistic),  so  that  N  =  Sc/Sw-  In  order  to  compute  the  inductance 
of  the  coil  we  use  the  formula  of  the  inductance  of  a  single  current  loop  multiplied  by  N  . 
Expressing  all  this  in  term  of  M,  R,  J  and  Sw  we  get : 
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where  RW,  MW  and  PW  are  respectively  the  wire  resistance,  mass  and  dissipated  power,  ?m 
is  the  wire  density  and  ?r  is  the  wire  resistivity.  A  sizing  could  be  the  following  : 
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In  the  presented  sizing  (M  =  I  6Am  ),  the  coil  radius  being  Im,  the  coil  will  probably  run 
around  the  spacecraft :  it  will  be  made  with  100  turns  of  6mm  copper  wire  for  a  total  length 
of  600m  and  a  weight  of  18Kg,  for  a  dissipated  power  of  34W. 


3.2.2  Tall  coil  with  ferromagnetic  core 

If  we  introduce  ferromagnetic  material  in  the  coil  core  the  sizing  is  slightly  different :  in  order 
to  be  sure  to  have  uniform  magnetisation  of  the  core  we  prefer  a  “tall”  coil.  This  increases  the 
mass  but  also  creates  a  stronger  field  in  the  coil  due  to  the  high  magnetic  permeability  of  the 
ferromagnetic  material  (pr).  Now  the  interesting  parameters  to  look  at  are  the  mass  and 
volume  of  the  magnetic  core  (Me,  Vc)  and  the  number  of  coils  per  meter  (n).  The  design 
parameters  are  M,  ?c  (core  density),  Bsat  (saturation  field  in  the  core),  Pr,  I  (the  current  in  the 
coils)  and  B  (the  ratio  between  the  coil  height  and  the  coil  diameter) : 


Vc  =  [i 


M 


B 


SAT 


*  SAT 
^SAT^ 


M 


^SAT 


M 

[Am^2 ] 

1, OOE+03 

Core 

vol . 

[m"3] 

6, 28E-04 

Core 

den . 

[Kg/in^3] 

8, OOE+03 

Core 

mass 

[kg] 

5, 03E+00 

Core 

Bsat 

[T] 

2, OOE+00 

n 

[1/m] 

3, 98E+01 

Mur 

[H/m] 

4, OOE+03 

L 

[H] 

5, OOE-03 

I 

[A] 

1, OOE+01 

R 

[m] 

2, 71E-02 

Beta 

(>2) 

[-] 

5, OOE+00 

H 

[m] 

2, 71E-01 

Wire 

res  . 

[Ohm.m] 

1, 70E-08 

Wire 

den . 

[Kg/m^3] 

8, 90E+03 

N 

[-] 

1, 08E+01 

Wire 

section  [m^2] 

3, OOE-06 

Lw 

[m] 

1, 84E+00 

Mw 

[Kg] 

4, 92E-02 

Rw  [Ohm] 

1, 04E-02 

Pw 

[W] 

1, 04E+00 

E  2 

In  this  case  (M  =  I  3Am  ),  the  magnetic  core  has  a  5Kg  mass  (27cm  height  and  2.7cm  radius) 
and  the  wire  dissipated  power  will  be  IW  :  with  respect  to  the  previous  case,  size  and  power 
consumption  are  much  smaller  (hence  we  can  hope  to  produce  larger  M  within  reasonable 
mass/power  limitations). 


3.3  Hub  magnetic  moment  sizing 

For  the  hub  magnetic  moment  we  present  here  two  solutions  based  on  superconducting  coils. 
In  order  to  obtain  the  rotating  field,  one  possibility  is  to  have  a  constant  magnetic  moment  on 
the  hub  and  spin  the  hub  itself.  Another  possibility  is  to  have  a  rotating  magnetic  moment,  but 
this  seems  complicate  to  obtain  with  superconducting  coils. 

3.3.1  Long  superconducting  coil 

Superconducting  materials  allow  large  currents  to  flow  without  any  resistance.  They  lose  their 
properties  above  critical  values  of  temperature,  current  density  and  magnetic  field  (the  three 
quantities  define  a  critical  surface).  We  consider  a  “tall”  empty  cylinder  (outer  radius  Ro, 
inner  radius  Ri).  The  magnetic  field  in  the  middle  of  the  cylinder  (minimum)  is  given  by  : 
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Expressing  everything  as  function  of  M,  Bq,  B  and  J  we  get : 
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In  this  sizing  (M  =  1  6Am  ),  the  cylinder  has  an  outer  radius  of  294mm,  an  inner  radius  of 
292mm  (radial  thickness  is  almost  2mm)  and  is  1.17m  long,  with  a  mass  of  30Kg. 


3.3.2  Flat  superconducting  coil 

In  the  case  of  a  “flat”  superconducting  coil,  the  sizing  is  done  just  as  a  classical  coil.  The  only 
difference  (on  top  of  the  fact  that  the  resistance  is  0)  is  that  the  field  in  the  core  cannot  be 
considered  uniform  (as  in  the  case  of  a  long  coil)  :  since  it  is  maximum  close  to  the  “wire”, 
we  need  to  check  that  it  value  does  not  exceed  the  critical  value  of  the  material. 


M 

[Am^2 ] 

1, OOE+06 

Sc 

[m-2] 

3, 18E-03 

J 

[A/m^2] 

1, OOE+08 

r 

[m] 

3, 18E-02 

R 

[m] 

1, OOE+00 

N 

[-] 

1, 06E+03 

Wire 

res  . 

[Ohmm] 

0, OOE+00 

Lw 

[m] 

6, 67E+03 

Wire 

den . 

[Kg/m^3] 

8, 90E+03 

Mw 

[Kg] 

1, 78E+02 

Wire 

section  [m^2] 

3, OOE-06 

Iw 

[A] 

3, OOE+02 

Rw 

[Ohm] 

0, OOE+00 

Pw 

[W] 

0, OOE+00 

Li 

[H] 

4, 79E-06 

L 

[H] 

5, 39E+00 

Mag  energy  [J] 

2, 43E+05 

Alpha 

[-] 

1, 07E+00 

Beta 

[-] 

6, 58E-02 

BO 

[T] 

5, 25E-01 

Bw 

[T] 

2, OOE+00 

E  2  2 

In  this  case  (M  =  1  6Am  )  the  Im  radius  coil  is  built  with  1000  turns  of  3mm 
superconducting  wire:  the  wire  is  6Km  long  and  weights  178Kg.  The  magnetic  field  in 
proximity  of  the  wire  is  2T. 


4  Conclusions 

The  proposed  actuation  system  is  challenging  :  it  requires  relative  position  and  attitude 
determination,  it  provides  “only”  average  force  and  couple,  it  delivers  instantaneous 
components  to  be  filtered  out  by  the  inertia.  It  is  also  very  challenging  on  the  technological 
side  :  having  large  magnetic  moment  on-board  may  force  the  introduction  of  constraint  in 
materials  and  on  the  integration  phase.  The  use  of  superconducting  materials  also  requires 
cryogenic  technology  and  has  many  unexplored  difficulties  (how  to  start/stop  the  large 
currents  circulating  in  the  coils,  safety  issues,  cost  and  total  mass  trade-offs).  Anyway, 
magnetic  actuation  could  be  an  alternative  to  other  micro  actuators  for  precision  formation 
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flying  mission:  their  main  disadvantage  is  the  1/r  dependency  of  the  magnetic  field,  limiting 
the  range  of  operation  in  terms  of  force/torque  and  distance.  But  compare  to  other  micro 
actuator  base  on  thrusters  they  offer  fully  scaleable  actuators  by  design,  they  don't  create 
accommodation  constraint  due  plume  effect,  contamination,  charging  effect,  etc . . . 
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